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■i ■ We present a brief account of the most salient properties of vortices in dilute 

r^ . atomic Fermi superfluids near a Feshbach resonance. 
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> ■ 1. INTRODUCTION 
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Since the creation of the first degenerate Fermi gas by DeMarco and Jin 

of K atoms and the optical trapping of ^Li atoms by Thomas' groujP'the 

^D ' last year and a half has produced an enormous experimental progress in the 



study of dilute atomic Fermi gases. The ability to manipulate the strength 
of the interaction by means of the Feshbach resonance opened extraordinary 

Cd I opportunities both from the experimental and theoretical points of view. 

There is no doubt in anybody's mind that dilute atomic Fermi clouds near 
1-^ . a Feshbach resonance should become superfluid, with a pairing gap of the 

Ch ! order of the Fermi energy at sufficiently low temperatures. The challenge 

O ■ from the experimental point of view is to realize these superfluids and espe- 

cially to demonstrate unambiguously the onset of superfluidity, or at least 
the formation of a condensate and be able to study its properties. A real 
/\ ' breakthrough was the creation and the subsequent study of the expansion 

j^ ■ of a strongly interacting degenerate Fermi gaP^. After that experimentalists 

have been able to study the formation of extremely weakly bound molecules 
(which we shall often refer to as dimers)^, the decay properties of ensembles 
of such dimera^, the BEC of dimers^, a number of features of the BCS to 
BEC crossover , the collective oscillationa^, the formation of some kind of 
condensate, with some still unclear propertiea^ and finally the appearance 
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of a gap in the excitation spectruirPl. 

Leggett and others have envisioned theoretically such a BCS to BEC 
crossovei!^ and were able to describe qualitatively its main features. Quali- 
tative features of the BCS dilute atomic Fermi superfluid have been discussed 
by a number of authors in recent yeara^^. The theoretical description was 
based essentially on the weak coupling BCS formalism, which is known to 
over predict the value of the gap by a significant factor^. The crossover the- 
ory of Leggett and its followers was based on a more or less straightforward 
extension of the weak coupling BCS formalism to the strong coupling regime. 
In the BEC limit there is an equally significant correction of this results^. 
As it was noted by BertsclP^, a dilute Fermi system acquires universal prop- 
erties at, what nowadays we call, the Feshbach resonance. The initial stud- 
ies of the Bertsch MBX challenge showed that such a system is stable^^'^. 
Only relatively recently that was confirmed both theoreticall} *^^ * ^^ -*^ and 
experimentallj'^. 

The discussion of some general properties of these systems, the charac- 
ter of the collective oscillations in trapped dilute atomic Fermi gases near a 
Feshbach resonance, which we planned to cover as well, along with the dis- 
cussion of the properties of atom-dimer mixtures, see Refs. I2()| are skipped 
here due to space limitations. 



2. Superfluid LDA and the Vortex State 

From the theoretical point of view the challenge is to be able to pre- 
dict and describe in a controllable manner the properties of these systems, 
which necessitates the development of accurate theoretical tools. We have 
shown recently how to extend the density functional theory ^^ to superfluid 
fermion systems, by creating the so called Superfluid LDA (SLDA)22I221, Jn 
the case of SLDA one needs to know the dependence of the energy density 
as a function of both normal and anomalous densities, unlike LDA when 
only the dependence on the normal density is sufficient. The rather accu- 
rate calculations of Refs. I17I18I allows us to construct the energy density 
functional (EDF) in the case of infinite homogeneous systems. Using the 
recent results of Chang et al. in Ref. E]we decided to re-parameterize this 
EDF and extend that parameterization away from the Feshbach resonance. 
In terms of the single quasi-particle wave functions, which define the nor- 
mal and anomalous densities the EDF of a superfluid system in the SLDA 
approach has the form: 
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(r)=^K(r)|2, r(r)=^|V^„(r)|2, u{v) =Y^vl{vW{v) 



where the spin degrees of freedom have been suppressed for the sake of sim- 
phcity. The dimensionless functions I3{x) and 7(2;) can be easily determined 
using the recent the results of Ref. il8i The SLDA equations have the same 
formal structure as the Bogoliubov-de Gennes equations for a system with 
density dependent contact interaction. In spite of the formal resemblance 
these equations describe such systems exactly, beyond the meanfield approx- 
imation, namely the ground state energy and number density distribution. 

The local anomalous density u{r) is unfortunately strictly speaking a 
diverging quantity and its evaluation requires a well defined regularization 
procedure. The evaluation of the pairing field also requires a renormal- 
ization procedure and so does actually the evaluation of the total ground 
state energy. The principles of the regularization and renormalization pro- 
cedures have been described by us previously, see Refs. I22|23l In the region 
of the Feshbach resonance, where the pairing gap is of order of the Fermi 
energy £p = ^^(37r^n)^'^/2?TT, these procedures have to be changed some- 
what, in order to ensure a better convergence of various quantities, number 
density, pairing field, total energy, etc. and these details will be presented 
elsewhere^. It suffices to add that unlike a number of other methods sug- 
gested recently in literature, aimed at dealing with similar divergences, see 
Refs. ISniand discussion in Ref. I26| and which require an active Hilbert space 
of a size 10^ . . . 10^ A^, in our approach the size of the active Hilbert space 
is typically much less than lOA^, where N is the total number of particles. 

The form of the EDF presented above is not unique to a certain extent. 
In particular one could have considered a density dependent mass in the 
kinetic energy term, as is often done in nuclear physics. Until full many- 
body calculations of the homogeneous matter of the kind described in Refs. 
I17|18|lf?l will provide more detailed information about the properties of such 
systems, beyond the energy per particle and the pairing gap as a function 
of the parameter \/na? , there is no unambiguous way one can determine 
whether the effective mass is different or not from the bare mass. There is 
no spin-orbit coupling as well and so far there is no compelling argument to 
expect its presence. However, one can expect gradient terms, in particular 
a dependence of the EDF on Vn(r). In principle such terms could be eval- 
uated, but their role is not expected to be ever dominant, though it could 
be significant. Phenomenologically, in nuclear physics it is established that 
such terms are quite important. The reason they are important in nuclear 
physics is because the radius of the interaction is comparable with the Fermi 
wave length. This is definitely not the case of dilute atomic gases for which 
nrfj <C 1 always (here tq is of the order of at most 100 A or so, the so 
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Fig. 1. The upper panel shows the number density profile around a vortex 
core, while the lower panel shows the profile of the pairing field. In the upper 
panel various curves correspond to l/kpa = (—0.5, —0.1, 0, 0.1, 0.3), from 
the highest to the lowest respectively. The order is reversed in the lower 
panel. We show also the (approximate) density profile of a vortex in a Bose 
dilute atomic gas of the same number density, with a dashed line for the case 



na" 



10 and with a dot-dashed line for the case na'^ 



10- 



called van der Waals length). One can come up with a similar qualitative 
argument in favor of an effective mass close in value to the bare mass, which 
was implicitly assumed by us. It is worth noting also that the present EDF 
is quite distinct from others suggested recently in literature, see for exam- 
ple Ref. 1271 which are typically based on one or another incarnation of the 
crossover model due to LeggetlP^. 

At this point all the elements needed, in order to perform a full self- 
consistent calculation of the vortex properties in a dilute atomic Fermi gas 
near a Feshbach resonance, are known. Without dwelling into technical 
detaila^^lSil -^g shall briefly discuss the most salient features. In Fig. 1 we 
show the number density profile n{p)/n{oo) and the profile of the pairing 
field A(p)/ei? around a vortex core, where p = {x,y) and the z-axis is along 
the vortex core. At the Feshbach resonance the asymptotic value of the 
pairing field is approximately one half the free Fermi energy of the free gas. 
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The actual Fermi energy at the Feshbach resonance is ~ O.Sep. As one can 
see from the lower panel in Fig. 1 the pairing field has a rather dramatic 
dependence on the scattering lengtlP^. 

What came as a big surprise before^^l and it is confirmed by the present 
results, based on a more accurate EDF, is the unexpected appearance of the 
prominent vortex core number density depletion shown in the upper panel 
of Fig. 1. For comparison we show there as well the number density profile 
of two vortices in dilute Bose gases of the same density and two different 
values of the corresponding atomic scattering length. As one can see the size 
of the vortex core is only 2-4 times smaller in the Fermi case than the vortex 
core in the BEC case. This fact let us conclude that a direct visualization of 
vortices should be easily achievable in the case of atomic Fermi superfluids 
around the Feshbach resonance. 



3. Conclusions 

The onset of superfluidity should undoubtedly be demonstrated by ex- 
citing and putting in evidence a superflow, and vortices are just about the 
only such modes in which superflow can be seen in such systems. Moreover, 
these vortices are also expected to form an Abrikosov lattice. We presented 
an analysis, based on a newly developed extension of LDA to superfiuid sys- 
tems, Superfiuid LDA (SLDA), of the properties of vortices at and near a 
Feshbach resonance. Surprisingly, like in the case of Bose superfiuids, vor- 
tices in Fermi superfiuids near a Feshbach resonance share common features, 
which should make them easily detectable. Vortices develop a pronounced 
density depletion in the core of a size comparable to the size of a core in a 
dilute atomic Bose gas. 
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